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PREFACE 


The  models  and  algorithms  in  this  paper  represent  a  synthesis  of 
results  from  different  works  by  several  different  individuals  over  the  past 
few  years.  In  particular,  I  would  like  to  acknowledge  the  contributions 
of  the  following  people:  Dr.  Walter  Chesnut  and  Dr.  Charles  Rino  of  SRI 
International,  Dr.  K.  C.  Yeh  and  his  colleagues  at  the  University  of  Illinois 
at  Urbana,  Dr.  David  Sachs  of  Science  Applications  Incorporated,  Dr.  Roy 
Hendrick,  Mr.  Robert  Bogusch,  Dr.  Dennis  Knepp  of  Mission  Research  Corporation, 
Dr.  Clifford  Prettie  of  Berkeley  Research  Associates,  and  Dr.  E.  J.  Fremouw 
of  Physical  Dynamics  Incorporated. 
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RADIO  WAVE  PROPAGATION  IN  STRUCTURED  IONIZATION 
FOR  SATELLITE  APPLICATIONS 


1.  INTRODUCTION 

The  fundamental  objectives  of  the  DMA  radio  wave  propagation  com¬ 
munity  have  been  to  determine  a  necessary  and  sufficient  representation  of 
structured  ionization  for  propagation  calculations,  to  characterize  the 
resulting  scintillated  signal  structures,  and  to  develop  methods  to  make 
available  scintillated  signals  for  systems  analysis  and  test  purposes. 

These  objectives  have  been  met  for  most  practical  applications. 

It  was  shown  in  References  1  and  2,  that  given  a  sufficiently 
thick  scattering  layer,  the  ionization  fluctuation  power  spectral  density 
was  adequate  for  most  propagation  purposes.  This  result  was  very  significant 
because  power  spectral  densities  can  be  developed  from  first  principles  or 
from  measurements.  If  a  more  complex  representation  had  been  required,  then 
it  would  have  been  doubtful  that  the  environment  problem  could  have  been 
adequately  solved  particularly  for  these  cases  where  empirical  characteri¬ 
zation  from  measurements  is  not  possible. 

The  characterization  of  typical  scintillated  signal  structures 
from  the  point  of  view  of  system  performance  evaluation  has  proven  particu¬ 
larly  simple  and  useful.  The  dominant  threat  to  systems  both  in  time  and 
space  is  from  Rayleigh  scintillated  signals  characterized  by  a  generalized 
power  spectrum  which  is,  in  turn,  parameterized  by  a  small  number  of  variables. 
This  simple  characterization  permits  a  separation  between  detailed  environ¬ 
ment  considerations  and  system  performance  evaluations.  Put  another  way. 


the  "environment"  as  seen  by  a  system  has  been  reduced  to  a  small  number  of 
variables  where  all  the  geometric  and  other  complexities  of  environments 
have  been  accounted  for.  Systems  evaluations  parameterized  by  the  reasonable 
ranges  of  these  signal  structure  parameters  thus  encompass  all  reasonably 
possible  environments  without  calculating  explicitly  all  possible  scenarios. 
The  appropriate  parameters  also  constitute  "environmental"  criteria  without 
the  necessity  of  dealing  with  specific  situations. 

The  generalized  power  spectrum,  because  of  its  properties,  permits 
easy  reconstruction  of  explicit  signals.  These  signals,  in  the  form  of 
computer  data,  may  be  used  to  exercise  software  models  of  satellite  systems. 
The  same  data  can  also  be  used  to  drive  test  hardware  for  testing  of  satellite 
hardwire  against  scintillated  signals.  The  remainder  of  this  paper  will 
develop  the  basic  models  and  algorithms  used  to  implement  the  concepts  just 
discussed. 


2.  ENVIRONMENT  CHARACTERIZATION 

The  first  step  in  any  propagation  study  is  to  characterize  the 
ionization  or  equivalently  the  index  of  refraction  fluctuation*  of  the 
ionospheric  propagation  environment.  Figure  1  illustrates  the  geometry  of  a 
typical  satellite  link.  The  index  of  refraction  fluctuations,  represented 
schematically  by  the  curved  lines  are  typically  elongated  along  the  magnetic 
field,  t  is  a  unit  vector  along  the  magnetic  field  ami  is  a  slowly  varying 
function  of  space  since  the  fie’d  lines  are  curved.  The  z  axis  and  the 
z  unit  vector  represent  the  propagation  line  of  sight.  The  transmitter 
by  definition,  is  at  2*0  and  the  receiver  is  at  t*  i^.  The  r  and 
s  unit  vectors  complete  with  t  an  orthogonal  coordinate  system  that  is 
used  to  define  the  structure.  The  orientation  of  r  and  s  is  chosen 
to  best  represent  any  anisotropy  of  the  index  of  refraction  structure 
about  the  field  line.  As  with  t,  r  and  s  may  be  slow  functions  of 
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position.  In  the  r,  s,  and  t  system  the  index  of  refraction  fluctuations 
are  often  represented  by  a  power  law  power  spectral  density. 


PSD(Kr,Ks,Kt) 


8iT3/2An?LrLsLtr(n)/r(n-3/2) 

(l*L2K2*L2K2*l2K2)n 


(1) 


where 

Kr,Ks,Kt  =  spatial  wave  numbers 

Lr,Ls,Lt  =  structure  outer  scales 
2 

An.  =  index  of  refraction  variance 
T(n)  3  gamma  function  of  argument  n 

2n  -  2  =  fluctuation  spectral  index 

i-i  > 

Both  in  situ  measurements  and  theoretical  considerations  imply  that 

n«2  for  ionospheric  structured  environments.  This  spectral  slope  is 

believed  to  be  caused  by  the  existence  of  very  sharp  edged  structures.  The 

fluctuation  power  spectrum  is  considered  a  locally  homogeneous  quantity  as 

10 

discussed  by  Taturski  .  This  means  that  the  parameters  can  be  slowly 
varying  functions  in  any  space  direction  with  respect  to  the  spatial  scale 
in  that  direction. 


The  structure  variation  of  the  index  of  refraction  perpendicular 
to  the  £  axis  dominates  the  propagation  effects  while  the  variation  parallel 
to  the  s  axis  enters  only  through  the  strength  of  the  integrated  phase 
variance.  Thus  Equation  1  must  be  transformed  to  a  frame  with  one  axis  being 
the  £  axis.  These  new  axes  are  defined  by 


m  a  A  a  a  . 

X  «  t*£/|t«2| 

(2a) 

(2b) 
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Two  rotations  suffice  to  accomplish  the  transformation.  First,  the  r  axis 
is  rotated  about  the  t  axis  by  angle  <p  to  become  parallel  to  the  x 
axis.  $  is  defined  by 

tsinifc  =  r  *  (t  *  z)  (3a) 

cos$  =  r  •  (t  *  2)  •  vb) 

Next,  the  t  axis  is  rotated  about  the  x  axis  by  angle  $  into  ttHe 
axis.  $  is  defined  by 

xsin$  =  t  *  2  (4a) 


eos<&  =  t  •  2  (4b) 

These  transformations  can  he  simplified  by  defining  new  effective  scale  sires 

L2  =  ijcos2*  ♦  l2sin2$  (Sa) 

L2  =  (L2sin^«J>  ♦  l2cos2$)cos2$  ♦  l2sin2$  (5b) 

y  r  s  t 

l2  o  (L~sin2$  ♦  Lfeos^Jsin'S  ♦  l^cos2$  (Sc) 

2  r  s  t 

Lxy  °  C Ls' LP cosCcosds inO  (5d) 

L  »  (l2-Lf)sin«cos$sin$  (Sc) 

L  s  (L2-L2$in2$-l*c0s*>$)cos$sin$  (5f) 

The  final  result  for  the  locally  homogeneous  index  of  refraction  power 
spectral  density  is 

S*5'2AnjL  t  L  r(n)/F(n-3/2) 

*  ■  o  ♦  •  s>y  *  w  * 2 v»v  * 2 w-.  * 21 w 


(6) 


N  > 


3. 


DERIVATION  OF  SIGNAL  STRUCTURE  PARAMETERS 


The  desired  signal  structure  parameters  are  derivable  from  the 

differential  phase  spectrum.  As  shown  in  Appendix  A,  this  spectrum  can  be 

2 

attained  by  multiplying  Equation  6  by  K  ,  where  K  is  the  radio  signal 
wave  number,  and  by  setting  Kz  equal  to  zero.  The  differential  phase 
spectrum  is 


dP.(K  ,K  )  8TT3/2K2An2L  L  L  r(n)/r(n-3/2) 
<b  x  yJ  i  r  s  t  v  J  K 


<t>  x’  y 
dz 


(1  +  L2K2  +  L2K2  +  2L  K  K  )" 
xx  y  y  xy  x  y 


The  total  mean  square  integrated  phase  fluctuation  is 


t  0°  c 

H  5/ 


dK  dP, (K  ,K  ) 
y  d> v  x5  yJ 

2ir  ~  dz 


o 


z 

■/- 


t  2 
no, 
<P 


where  do'/dz  is 
<P 


%  _  2TT1/2r(n-l)  2,  2  m 

~Sz  '  TTn-'372j  K  i  z  {  ) 

7  2  2  \ 

l  =  L  L  L./(L  L  -L  )  (10) 

z  r  s  t  v  x  y  xyJ 

l_  is  the  effective  structure  scale  along  the  z  axis.  The  spectrum  in 
z 

Equation  7  can  be  rewritten  as 


dP.(K  ,K  ) 
d>  x  y 


do2  4ti(L2L2-L2  )1/2(n-l) 

_ |  v  x  y  xy  _ 

dZ  (1  +  L2K2  ♦  lV  *21  Kli  )n 

v  xx  y  y  xy  x  yJ 


xy  x  y 
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Equation  11  can  be  Fourier  transformed  to 


<V*,y)  <*>2  pb'Vi“» 

<“  *  <2  2n"2r(n-l)  '  ’ 

where  7  -  ~  - 

2  y-yy 

L2L2  -  L2 
x  y  xy 

Equations  9,  11,  and  12  are  the  fundamental  equations  from  which  all  of 
the  required  propagation  quantities  can  be  derived. 

The  first  propagation  quantity  is  the  mean  square  log  amplitude 

'  a 

fluctu;:  -o.:  This  quantity  is  used  as  a  switch  to  determine  the  onset  of 

Rayleigh  scintillated  signals. 


-*)*  diyw 


o  —CO  -Cd  .  . 

(d.|.2)(2  -!)I 

Mf:)5  — Wt-  '•  : 

■  '  '  V 

The  mean  square  log  amplitude  fluctuation  can  be  approximately  expressed  by 
(see  Appendix  S) 


■2= iHrs^rV 

X  ~  2  J  Us  d‘ 


Equation  14  is  accurate  to  about  twenty  five  percent  everywhere  and  usually 
to  within  ten  percent  when  M{*}  is  always  small .  The  applicable.  range 
of  n  is  1.7  s  n  5  2.75. 
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The  remaining  signal  structure  parameters  are  derivable  from  the 
mutual  coherence  function  of  U(z,r,K)  where  the  propagating  carrier  signal 
is  expressed  as 


2  _2, 1/2 


Ufz  ~  KleiK(ir|*"  +  z“:) 
E(z>r>K)  =  U(Z^,K^e - - - 


The  mutual  coherence  function  for  U  is 


(15) 


G(z,x,y,n)  U*(z,r1 , K^UCz ,r2«K2) 


(16) 


where 

n  =  (K2-K1)/K2 
xx  +  yy  =  r2  -  r  x 
The  equation  for  G  is 


where 


(17) 


The  differential  total  integrated  phase  variance,  do^/dz  is  evaluated  at 
K^Kj.  Equation  17  includes  the  effects  of  spherical  waves  which  permits 
the  transmitter  and  the  receiver  to  be  located  arbitrarily  with  respect  to 
the  scattaring  medium. 


Let  n  equal  zero. 


Then  Equation  17  can  be  written 


(18) 
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with  the  formal  solution 


where 


G(z,x,y)  =  e 


G(0,x,y)  =  1  . 


¥•?)]* 


The  bracket  in  Equation  19  can  be  approximated  for  »  1  by 

s' 


where 


m  =  2 ,  n  s  2 


m  =  2n  -  2,  n<2 

B  =  constant  of  order  unity 


(19) 


(20) 


The  signal  decorrelation  distance  is  defined  as  the  e  1  point  of 

G.  Thus,  for  o'?  »  1 , 

$ 

zt 

r\  2  m 

f  if  B n  {if  )  4*  “  1  <21) 

o 

defines  p^,  the  dccorrolati«\.  listance.  The  decorrelation  distance  is  a 
function  of  direction  in  the  x-y  plane.  The  signal  decorrelation,  time  is 
defined  as  the  decorrelation  distance  in  some  direction  dividcv'  by  the 
effective  velocity  in  that  same  direction.  The  effective  velocity  is 


where 


=  transmitter  velocity 


V  =  index  of  refraction  structure  velocity 


Let 


Vye  =  receiver  velocity 


From  Equations  20  and  21,  the  perpendicular  decorrelation  time,  ta,  is 


where 


(23) 


B(n.O^)  =  minimum  ,  (-0.  J4n2  ♦  2. Sin  -  2.00)  (o*Bn) 1/ra  j 

7 

The  B(n,o”)  function  is  used  to  insure  the  proper  limit  on  for  small 
o2.  The  decorrelation  time,  Tj ,  is  with  respect  to  the  x-y  plane.  The 

'  A 

signal  also  decorrelutes  with  time  for  motion  along  the  z  axis  as  a 
result  of  the  angular  scatter  of  the  energy.  We  will  assume  for  the 
remainder  of  this  paper  that  n  2  and  o“  »  1  to  continue  the  development. 
Restricting  n  limits  the  classes  of  environments  but  allows  a  more  general 
geometry . 


Let  us  define  a  set  of  fixed  axes,  u  and  v,  perpendicular  to 

A  AAA  a 

the  z  axis  such  that  u*v»j.  At  each  point  on  the  z  axis,  the  rotation 
angle  between  the  x  and  the  u  axis  is  defined  by 

zsin9  =  x  *  u 

A  A 

COsO  a  X  •  U 

14 


Let 


(29a) 


n 


C  a  ~  |c  +c  +((C  -C  )2+4C2  )1/,21 
p  2  [  u  v  u  u  V  uv'  J 

C  =  \  fc  +C  -((C  -C  )2+4C2  )1/21 
q  2  I  u  v  1  *■  u  v'  uv'  j 


Then 


-<V2*V2) 


G(zt,P,q)  =  e 


The  parallel  signal  decorrelation  time  is  defined  as 


3.7K 

(C'  v')  *,^2T  (V  -V  J  •  z 

p  q  '  fV  re  st' 


(29b) 


(30) 


T“  *  -.2/3^2/3,3/21^  T? 


(31) 


The  final  decorrelation  time  is 


tq  =  minimum^,  ty) 


(32) 


The  signal  power  spectrum  that  is  most  consistent  with  the  proceeding 
development  is 

2 

(33) 


-(TrfT  y 

P(f)^re  ° 


S  '  o 

Occasionally,  the  following  form  is  used  to  be  conservative  when  n  =  2 


vt  /l .66 

Km  •  — 2 — i\vn 


TeT 


(34) 
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Equation  25  can  be  Fourier  transformed  resulting  in  the  energy  angle  of 
arrival  spectrum.  -  .. 

"  **-•"*  a  9  r 

*  U  V  V  ■ 


VW  a 


d**2c  a  a  »c 

K  uv  uv  u  v  v  u 
, ‘  4  I  CC-C2  J 

L  U  V  UV  J 


4ff(C  C  -C*  ) 
v  u  V  uv7 


T7T 


(3S) 


I 


I 

.r-- 


t 

f 

% 


or  in  the  p-q  frame 


WV 


~;v  ft 


4*(CPV 


1/2 


where  both  angular  distributions  are  normalized  such  that 


/d9v/ d9u  WV  ■  1 
-00  00 


00  00 

/de  fde  pa(8  ,e )  =  i 
p  J  q  0V  p  q 

-co  oo 


(36) 


(37a) 


(37b) 


The  final  parameter  to  be  calculated  is  the  variance  of  the  signal 
energy  arrival  time.  Using  a  formalism  originally  developed  by  Yeh  and  Liu 
the  effective  energy  arrival  time  variance  is  defined  as 


where 


o  o 
z r 
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o^p  is  the  Rayleigh  phase  variance  derived  in  Appendix  B.  It  is  used  to 
insure  that  0^  includes  only  the  time  delay  effects  associated  with  the 
Rayleigh  portion  of  the  signal.  The  frequency  selective  bandwidth  is 
defined 


f0  =  l/2not  (39) 

This  completes  the  development  of  the  basic  signal  structure  parameters. 
In  the  following,  these  parameters  will  be  used  to  define  the  generalized 

i. 

signal  power  spectrum. 

4.  THE  GENERALIZED  POWER  SPECTRUM 


Let  us  assume  that  the  thickness  of  the  scattering  layer,  zs>  is 
small  compared  to  the  distance  from  the  layer  to  the  receiver  and  to  the 
transmitter.  This  is  modeled  mathematically  by  letting 


4<S(z"2i) 


(40) 


where  6(z)  is  the  pirac  delta  function  and  z^  marks  the  center  of  the 
scattering  layer.  Then 


t  The  use  of  this  expression  throughout  constitutes  the  Delta  layer 
approximation  which  is  the  basis  of  a  set  of  simple  propagation 
expressions  for  well  defined  scattering  regions  (sec  Appendix  C). 
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£p,  f^,  and  have  implicit  dependence  on  f2>  In  actual  appl ication'- 

the  system  bandwidth  is  usually  small  compared  to  the  carrier  frequency,  fc 
Thus  the  above  parameters  vary  iittle  over  the  bandwidth  and  as  it 
appears  everywhere  in  Equation  31,  except  in  Af,  can  be  replaced  by  f..  This 
is  consistent  with  dropping  the  1  -  n  terms  earlier.  With  these  assumptions. 
Equation  39  describes  a  statistical  process  that  is  stationary  ir.  both  space 
and  frequency. 


The  delta  layer  approximation  is  not  always  justified  particularly 
when  the  scattering  layer  is  thick  and  comparable  with  the  other  problem 
dimensions.  Appendix  0  demonstrates  that  finite  layer  thickness  can  be 
compensated  for  by  correcting  fp  and  P  with  a  multiplication  factor 
and  by  multiplying  G(zt,p,q,Af)  by  exp(iR’Af)  whore  R*  is  a  constant 
independent  of  Af.  This  term  can  be  neglected  because  it  results  in  a 
uniform  time  shift  of  the  signal.  These  corrections  apply  for  all  cases 
ranging  from  the  thin  layer  to  the  case  where  both  the  transmitter  and 
receiver  are  buried  in  the  structured  layer.  The  success  of  the  simple 
scaling  of  f^  and  P  demonstrates  the  generality  of  the  form  shown  in 
Equation  41  assuming  fp  and  f^  are  properly  chosen. 
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Equation  41  can  be  specialized  to  two  cases  which  bound  the 

frequency  selective  effects.  The  first  case,  where  Cp  =  C^,  assumes 

isotropy  of  the  effects  about  the  propagation  line  of  sight.  In  this  case 

only  one  coordinate  can  be  relevant  so  q  can  be  set  to  zero.  Finally, 

2  2  2  2 

the  remaining  spatial  argument,  p  transforms  to  At  fx^  when  the 
effective  motion  is  considered.  The  result  is 


G2(zt,At,Af)  = 


a2  2  2 

<*R/Afr  At i  l  1  \ 

;  2  ifcJ '  t2  \ 


where 


f if  uw*Ivw*<^i,'1/2 


2 

o^p  has  been  used  to  limit  G2(*t,At,Af)  to  the  Rayleigh  portion  of  the 
signal.  The  redefinition  of  f£  preserves  which  was  defined  in 
Equation  38. 


The  second  case,  whore  ^>>0^,  assumes  that  there  is  signal 
variation  in  only  one  direction  with  respect  to  the  propagation  line  of 
sight.  Transforming  the  spatial  dependence  to  time,  we  get 


GjCz^At.Af)  « 


(1-iAf/fp 


[-2  •  <j  » 

$R/Afr  At"  /  1  \ 

•  i  J  ■  ti  (rw?ri 
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where  fj  =  f£//2.  The  generalized  power  spectrum  is  the  two  dimensional 
Fourier  transform  of  G(zt>At,Af). 

00  00 

r(f'T)  =  /d(Af)  /  d(At)G(zt>At,Af)e"i2lrAfTei2irfAt  (44) 

where 

f  =  frequency  coordinate 
x  =  time  delay  coordinate 


Evaluating  Equation  44  with  Equations  42  and  43 


r2(f,T) 


where 


(45) 


*1 


(f.T)  -  2 


-<*V> 
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2 
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?K«T0f)2 


(46) 


Both  Tj  and  1^  arc  nondieensional  and  normalized  to  unity  when  integrated 
over  all  f  and  t.  When  integrated  over  all  t,  they  reduce  to  Equation  33, 
the  fiat  fading  spectrum.  Both  functions  depend  on  three  nondimens ional 
quantities:  TQf,  f’r,  and  ^/(f'o^)  which  must  be  specified  to  define  the 
generalized  power  spectra. 
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The  generalized  power  spectra  in  Equations  45  and  46  represent 
a  stationary  process  with  a  Rayleigh  distributed  amplitude.  This  consti¬ 
tutes  a  complete  statistical  description  of  the  signal  structure  excluding 
mean  and  asymmetric  bandwidth  effects.  Figures  2  through  8  show  selected 
plots  of  the  generalized  power  spectra.  For  o,Df'/f  large,  the  gaussian 
term  in  frequency  dominates.  For  o^f'/fc  small,  the  energy  arrives  at 
higher  frequencies  for  longer  delays.  This  behavior  is  intuitively  pleasing 
because  increasing  frequency  corresponds  to  increasing  energy  angle  of 
arrival  with  respect  to  the  propagation  line  of  sight  which,  in  turn, 
implies  longer  delay  because  of  longer  geometric  path  lengths.  For  (^(f.T), 
some  energy  is  at  low  frequency  regardless  of  the  delay  because  of  energy 
scattering  from  regions  out  of  the  plane  determined  by  the  line  of  sight  and 
the  direction  of  the  effective  velocity.  The  most  significant  difference 
between  r^f.T)  and  ^(f.x)  is  in  the  total  energy  arriving  at  any  delay 
as  shown  in  Figures  9-12.  T^ff.t)  is  usually  the  more  serious  threat  because 
it  leads  to  more  delay. 


The  effects  of  a  disturbed  linear  channel  can  be  expressed  in  terms 
of  the  channel  impulse  response  function,  h(t,t).  The  received  signal  modu¬ 
lation  is 


Sr(t)=  /  h(t.t)St(t-T)dT  (47) 

whore  S^Ct)  is  the  transmitted  modulation.  The  impulse  response  function  is 
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Figure  3.  Two-dimensional  generalized  power  spectrum. 
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Figure  7.  One-dimensional  generalized  power  spectrum. 


arrival  as  a  function  of  delay 


Figure  10.  Energy  arrival  as  a  function  of  delay. 


where  iu  is  the  mean  index  of  refraction.  U(zt>t,f)  is  the  carrier  signal 
defined  in  Equation  15  but  expressed  as  a  function  of  time  as  perceived  at 
the  receiver.  The  mean  index  of  refraction  can  be  expanded  about  the  carrier 
frequency  resulting  in  the  following  form  for  the  channel  impulse  response 
function. 
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(49) 

where  TEC  is  the  total  electron  content,  r  is  the  classical  electron 

2  0 

radius,  and  only  terms  through  (f-f£)  have  been  retained.  Equation  49  is 
a  rather  complete  description  including  doppler  effects  from  changes  in  z+ 
and  the  most  important  TEC  effects:  phase  shifting,  group  delay,  and  disper¬ 
sion.  Let  the  scintillation  impulse  response  function  be  defined  by 


h#(t.T) 


J df  U(zt,t,f)e'l2lTfT 


(50) 


Since  UCz^.t.f)  is  Rayleigh  distributed,  then  so  is  h  (t,x).  With  some 
manipulation,  it  can  he  shown  that  the  frequency  power  spectrum  of  hs(t,T) 
is  either  Equation  45  or  46.  Thus  Equation  45  or  46  provides  a  complete 
statistical  description  of  the  scintillation  impulse  response  function.  The 
mean  value  of  t  in  Equation  50  is 


or 


(51a) 
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(51b) 


depending  on  which  of  the  two  generalized  power  spectra  is  used. 
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Equations  45  or  46,  49,  and  50  permit  the  sampling  of  the  channel 
impulse  response  function.  First,  using  methods  described  in  Appendix  E, 
samples  of  hs(t,x)  are  generated.  For  the  general  case,  hg(t,x)  is  next 
Fourier  transformed  with  respect  to  T,  multiplied  by  the  dispersion  and 
time  delay  terms,  inverse  transformed  with  respect  to  f,  and  finally  multi¬ 
plied  by  the  doppler  and  TEC  phase  term.  In  most  cases  scintillation  effects 
are  evaluated  separately  from  the  other  effects  and  the  scintillation  impulse 
response  function  is  used  as  the  channel  impulse  response  function.  If  the 
scintillation  delay  is  small  with  respect  to  the  modulation  symbol,  then 
the  delay  can  be  ignored.  Then  Equations  45  and  46  can  be  replaced  by 
Equations  33  or  34.  The  scintillation  impulse  response  function  is  only  a 
function  of  time.  Thus, 


2irf  icr  TEC 


This  completes  the  discussion  of  the  generalized  power  spectrum 
and  its  value  in  providing  a  complete  statistical  description  of  a  scintil¬ 
lated  signal.  Most  important,  methods  were  presented  to  economically 
generate  disturbed  signals  for  test  and  evaluation  applications. 

5.  SUMMARY 

The  preceding  sections  have  developed  a  sequence  of  models  and 
algorithms  that  specify  the  index  of  refraction  fluctuation  environments, 
calculate  various  signal  structure  statistical  parameters,  and  lastly, 
permit  the  sampling  of  the  channel  impulse  response  function  consistent  with 
the  signal  structure  statistical  parameters.  The  geometric  aspects  of  the 
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problem  were  treated  in  all  generality  using  a  locally  homogeneous  approxi¬ 
mation.  All  of  the  preceding  development  applies  to  2sn<3.  This  range 

includes  the  presently  held  value  for  strongly  structured  ionospheric 

2 

environments.  We  usually  assumed  o^»lt  but  this  restriction  is  not  very 
serious.  It  can  be  neglected  without  seriously  impacting  any  results  as 
they  may  impact  communication  or  radar  systems. 
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APPENDIX  A 

DERIVATION  OF  THE  LINEAR  INTEGRATED  PHASE  SPECTRUM 


Let  the  integrated  phase  be  written  as 
4>(r)  =  K  [ai\.  (r,z)dz 


(A-l) 


where 


z  =  coordinate  along  the  propagation  line  of  sight  (LOS) 
r  =  coordinate  vector  in  perpendicular  to  the  LOS 


An^  =  index  of  refraction  fluctuation 


K  =  wave  number  of  carrier  signal 
Now  look  at  a  portion  of  the  integral  in  Equation  A-l. 

rzl^z 

A<j>(r, 


^'l , z  1 )  =  K J  &ni(r1,z')dz 


'1 


The  correlation  function  for  A4>  is 


z^+Az  Z2+Az 


%(p,£)  =  r/  dz 


Az 


=  rAz J  dnR(p,n-e)(i-!n|/Az) 
-Az 
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(A-2) 


j  f  ^Z?  ^ni^l,2^^ni^r2,z'^  (A-3) 


(A-4 ) 


where 


P  =  rl  -  r2 


e  =  z1  -  z2 


R(p,z-z')  =  An.(r1,z)Ani(r2>z’) 

If  we  assume  that  R(p,z-z*)  has  a  decorrelation  distance  L  for  the  second 

z 

argument,  then  it  is  clear  that,  as  Az  becomes  larger  than  Lz’  %  is 
significant  only  for  small  e.  As  Az  increases,  the  contributions  to  the 
total  phase  correlation  function  from  each  layer  of  thickness  Az  become 
uncorrelated.  Thus  for  Az  larger  than  and  assuming  that  the  statis¬ 
tical  parameters  in  R  do  not  vary  greatly  over  Az,  we  can  write  the  phase 
correlation  function  as 


00 

yp)  -  K2Jdz  J  de  R^Cp.e) 


(A-5) 


where  R^  depends  on  z  through  its  statistical  parameters.  Equation 
A-5  allows  us  to  define  the  differential  linear  phase  spectrum.  Let 


(Kp)  =  J d2p  elKp  PJ dz  jdt  K2RA<J((p,e) 


CA-6) 


Since  we  have  assumed  that  R^CP.e)  changes  slowly  over  Az  and  hence  Lz, 
the  integral  over  z  can  be  taken  to  the  outside  and  differentiated.  The 
result  is  the  differential  linear  phase  spectrum. 


The  operations  in  A-7  are  equivalent  to  multiplying  the  in  situ  index  of 
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refraction  spectrum  by  K  and  setting  the  wavenumber  along  the  LOS  equal  to 
zero.  Thus  for  a  generalized  power  law  spectrum 


dP»(ICX*V  8TT3/2K2An^LxLyLzr(n) 

dz  r(n-3/2) (l+L2K2+LyKy)n 


(A-8) 


In  the  derivation  of  Equation  A-8  it  was  assumed  that  Lz  was 

smaller  than  Az  and  hence  the  scattering  layer  thickness.  If  L  becomes 

z 

comparable  to  or  larger  than  zg,  the  layer  thickness,  then  the  earlier 
derivation  fails.  Let 


4>(r)  =  (i)  KJ  Cr,z) 


-  -fz/zs) 

r,z)e  dz 


(A-9) 


where  the  exponential  represents  a  finite  scattering  layer  along  the  LOS. 
Ani  is  for  present  purposes  taken  as  a  homogeneous  stationary  process. 
After  some  manipulation 


R,(p)  =  K2Zs  Jr(P, 


-e2/2z2 
e)e  s  de 


(A-10) 


where  global  homogeneity  has  been  assumed.  The  phase  power  spectrum  is 
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If  z  /L  is  large,  then  the  exponential  integrates  out. 
s  z 


P(K  ,K)  * 


8iT3/2An?K2L  L  L  z  T(n) 
i  x  y  z  s 


VV  Hxv2T 2  J  2.n 


(UK‘L%KyLp1‘r(n-3/2) 


(A- 12) 


Equation  A-12  is  Equation  A-8  integrated  over  z.  Now  if  z  / L  is  small, 

22  22  22  sz 

the  spectrum  has  two  regimes.  For  KL  +KL  >L/z  ,  Equation  A-12  still 

xx  yy  zs  22  22  22 

holds  since  the  exponential  dominates  the  integral.  For  K  L  +K  L  <L  /z  , 

xx  y  y  z  s 

the  power  law  dependence  on  K£  dominates  and 
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(A- 13) 


Thus,  the  actual  phase  spectrum  has  two  distinct  parts.  The  low  frequency 
part  is  one  power  shallower  in  slope  and  goes  as  the  square  of  the  layer 
thickness. 


In  general,  the  spectrum  in  A-12  is  most  accurate  as,  with  one 
exception,  Lz  doesn't  get  appreciably  larger  than  z .  For  these  cases 
the  use  of  A- 13  would  be  a  gross  overestimate  of  the  effects.  The  one 
possible  exception  is  when  the  LOS  is  nearly  parallel  to  the  magnetic  field 
in  the  ionospheric  scintillation  case.  In  this  case  the  outer  scale  along 
the  magnetic  field  may  be  very  long  compared  to  the  thickness  of  the  region 
with  any  scattering  strength.  Even  in  this  case,  Equation  A-12  will  not 
underestimate  the  effects  and  gives  the  best  representation  for  tho  larger 
transverse  wave  numbers. 
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significantly  different  in  size,  then  the  integral  is  proportional  to 

M(z)~^2.  Thus  the  second  term  on  the  right  in  Equation  B-4  does  not  con- 
2 

tribute  to  y  when  M(z)  is  large.  For  M(z)  large 


(B-5) 


2 

For  M(z)  small,  the  x  integral  is  dominated  by  values  of  K 


and  Ky  where  L2K2+L2Ky+2LxyKxKy  » 1 .  Then 
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First,  we  rotate  in  K  -K  space  to  eliminate  the  K  K  term. 
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where 
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The  first  bracket  is  approximately  equal  to  (n-l)/8.  The  last  bracket  is 
approximately 
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The  final  expression  for  x  with  small  M(z)  is 
~2  i  j  do2 

X  /  dzMC2)n_1g(n)  -4 


(B-9) 


(B- 10) 


where  g(n)  equals  0.68,  0.79,  and  0.83  for  valuos  of  n  of  1.75,  2.0,  and 
2. S  respectively.  For  simplicity  we  will  let  g(n)  be  unity.  The  error 
incurred  is  small  particularly  for  n  equal  or  greater  than  two  which  is 
the  most  common  case. 


result. 


This  last  case  can  be  combined  with  Equation  B-5  to  get  the  final 


mm 


Comparison  with  explicit  propagation  calculations  has  shown  Equation  B-ll 
to  be  accurate  to  about  twenty  five  percent  everywhere  and  within  ten  percent 
if  M(z)  is  small  over  the  entire  integral. 


For  small  M(z) ,  the  phase  variance  is  divisible  into  the  Rayleigh 
phase  variance  and  the  remainder.  The  Rayleigh  phase  variance  is  defined  by 


,2  a  & 


~  4n 


ill  f  da  f  dz 

-2j  a2n  lJ  2 


t  do2 


n-1 


(B-12) 


where  the  notation  of  Equation  B-8  has  been  used.  ac  is  chosen  by  using 

Equation  B-8  to  determine  that  portion  of  the  integral  over  a  that 

~2 

contributes  negligibly  to  x  • 


(n- 


-1)  f  sin2(a2)  j  f  dz  ^°<J>  w/  ^n-1 

Zz  J  a2n-l  i!1J  T~ST  Mw  -X 


t  .2 
do. 


4n 


o  o 

Equation  B-13  can  be  simplified  using  Equation  B-ll  and  setting 


sin2 (a2)  »  a4 


Integrating  over  a,  remembering  that  M(z)  is  small 


6-2n 


6-2n 


(ndV 


4n-2 


1  -l 


Similarly  with  Equation  B-12 


(2n 


J _ I Vjl  x2l 

-2}a2n-2  [4«“2  J 
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(B-13) 


(B-14) 


(B-1S) 


APPENDIX  C 

THE  DELTA  LAYER  APPROXIMATION 


This  approximation  consists  of  using 


da 


4>  . 


dz 


a^(z-Zl) 


(C-l) 


where 


,1/2 


_2  f\  2tt  '  T(n-l)  J2,  2  J 

a4> =  J  Lz  ~f TiT-W  K  toi dz 


Lz  = 


LrLsLt 


(LxLy’Lxy)1/2 


evaluated  at 


z 


and  Zj  marks  the  structured  layer.  In  the  following  all  of  the  scale  sizes 
used  are  evaluated  at  . 


where 


The  mean  square  log  amplitude  fluctuation  is 

t  ‘ir^r1 

x  =T[TTm^J 


M(Zj)  - 


‘f  = 


K(LxLy-L^ 


(C-2) 
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The  variance  of  the  energy  arrival  time  is 


4„2  2  ,A  .4 


,  0*  2a7B  zf  (L  +L  +2L  ) 

2  _  4>R  »  n  f  v  x  y  xyJ 

t_K2c2  K4c2  (L2L2-L2)2 

1  x  y  xyJ 

and  the  frequency  selective  bandwidth  is 
f0  -  (2™,)-1 


(C-7) 


(C-8) 


Finally 


f'  = 
2 


“ 

2  1 

_F_ 

/2 

f2 

L  o 

- 1 

— 

F 

i 

°<J)R 

A 

f2 

c  J 

1-1/2 


1-1/2 


(C-9) 


(C-10) 


where 


F  =  1  ♦  0.038R"2’4 

F  ■  l  ♦  0.14R’05 

R  =  ze/z 
f  s 


RsO.S 
R  >0.5 


2 

are  the  correct  terns  to  be  used  in  the  generalized  power  spectra.  a.u 

yK 

the  Rayleigh  phase  variance,  was  defined  in  Appendix  B.  The  correction 
factor,  F,  was  derived  in  Appendix  D  to  correct  for  the  finite  structured 
layer  thiekness, 


«»  * 
$ 
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The  various  propagation  cases  can  be  parameterized  by  any  two  of  the  following 
three  numbers. 


R1  ■  Vzs 


R2  "  (VZ1)/2S 

R 


Vzs 


R1R2/(R1+R2) 


Zj  is  the  coordinate  of  the  center  of  the  layer.  Equations  D-2  and  D-3  were 
solved  by  a  combination  of  analytic  and  numerical  techniques  for  the  values 
of  Rj  and  R2  shown  in  Figure  D-l.  The  results  show  that  G(p,Af)  can  be 
well  represented  by 

2 


Gfp.Af)  « 


exp(iGAf/f ') 


(1-iAf/Ff^) 


exp 


°WAf'2 


2\f 


)  -  6)  M 


iAf/Ffl 


CD-4) 


where 


p  = 
2 


2  2 
KcL 

X 


4itcr(j)zfBn 


2.4 


F  =  1  +  0.038R 

,-0.5 


F  =  1  +  0.14R 


R£0.5 
R  >  0. 5 


Equation  L-4  is  identical  with  Equation  42  in  the  text  except  for  the  F 
and  G  corrections.  The  G  term  can  be  ignored  as  it  represents  a  uniform 
time  shift. 


Figures  D-2  through  D-ll  compare  the  solutions  for  G(p,Af)  (symbols) 
with  Equation  D-4  (solid  curves).  The  excellent  agreement  is  evident. 

Figure  D-12  plots  the  correction  factor,  F,  as  a  function  of  R  that  best 
corrects  Equation  D-4. 


o-i  Sreenivasiah,  I.,  et  al.,  "Two-Frequency  Mutual  Coherence  Function  and 
Pulse  Propagation  in  a  Random  Medium:  An  Analytic  Solution  to  the  Plane 
Wave  Case,"  Radio  Science,  Vol.  11,  Number  1_0,  775-778,  October  1978. 
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Figure  0-3.  Comparison  of  mutual  coherence  functions 


functions. 


Figure  D-6.  Comparison  of  mutual  coherence  functions 


Comparl 


APPENDIX  E 

GENERATION  OF  EXPLICIT  SAMPLES  OF  RAYLEIGH  SCINTILLATED  SIGNALS 


The  statistics  of  the  signal  are  represented  by  the  generalized 
power  spectrum,  r(f,x),  which  is  an  implicit  function  of  f',  o^R,  and  T0« 
Since  the  process  represented  by  r(f,i)  is  stationary,  r(f,x)  can  be  thought 
of  as  a  colored  noise  spectrum.  This  means  that  samples  of  the  signals  can 
be  generated  parametric  in  delay,  x.  Let  a  signal  sample,  for  x^  delay,  be 


U(t,T.) 


JL 

2T 


N 

£ 

n=-N 


C  ft.) 
nv  i' 


.  mrt 


(E-l) 


where  2T  is  the  length  of  the  sample  and  the  "ange  of  t  is  -Tst^T. 
2N  is  the  number  of  discrete  frequency  points.  Lot 


T  >  100  t 

o 

T/N  <  tq/10 

The  Fourier  coefficients  in  U-l  can  then  be  assumed  zero  mean,  independent, 
and  gaussian  distributed  variables  with  a  variance  of 

x.nVr/2  (n*l/2)/2T 

■  (if/  *  / ,f  <e-2> 

X,-&t/2  (n«l/2)/2T 


where  ot 


is  the  delay  resolution  element. 


Each  Fourier  coefficient  can  be  sampled  by 


C  (t.)  =  C*(T.)C  (t.) 
iJ  n  1  i' 


1/2  ,  ,,  i2irR 

l-MRj)]1'  e 


(E-3) 


where  R^  and  are  two  random  numbers  from  a  distribution  whose  elements 
are  uniformly  distributed  between  zero  and  one.  The  sample  signal  of  delay 
is  found  by  evaluating  Equation  E-l.  This  process  is  repeated  for  all 
T. .  For  many  applications  this  may  be  the  final  step  because  U(t,T^)  is 
also  a  sample  of  the  channel  impulse  response  function. 


A  special  sampling  case  occurs  in  flat  fading  where  the  signal 
delay  is  small  compared  to  the  symbol  period.  The  argument  in 
Equation  E-l  can  be  dropped.  The  variance  of  the  coefficents  are 

_  (n+l/2)/2T  » 

C*Cn  =  (2T)2J  df  J  dTT(f,T)  (E-4) 

(n-l/2)/2T  -• 

The  coefficients  are  sampled  by  Equation  E-3. 
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